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Abstract. We establish the presence of a spectral gap near the real axis for the damped 
wave equation on a manifold with negative curvature. This results holds under a dy- 
namical condition expressed by the negativity of a topological pressure with respect to 
the geodesic flow. As an application, we show an exponential decay of the energy for all 
initial data sufficiently regular. This decay is governed by the imaginary part of a finite 
number of eigenvalues close to the real axis. 



1. Introduction 

One of the standard questions in geometric control theory concerns the so-called sta- 
bilization problem: given a dissipative wave equation on a manifold, one is interested in 
the behaviour of the solutions and their energies for long times. The answers that can be 
given to this problem are closely related to the underlying manifold and the geometry of 
the control (or damping) region. 

In this paper, we shall study these questions in the particular case of the damped wave 
equation on a compact Riemannian manifold (M, g) with negative curvature and dimension 
d > 2. For simplicity, we will assume that M has no boundary. If a G C°°(M) is a real 
valued function on M, this equation reads 

(1.1) {dt - Ag + 2a{x)dt)u ^ , {t,x)eM.xM, 
with initial conditions 

m(0, x) = ujq{x) e 
idtuiO,x) = LUi{x)eH°. 

Here H"^ = H^{M) are the usual Sobolev spaces on M. The Laplace-Beltrami operator 
Ag = A is expressed in local coordinates by 

(1.2) Ag = ^dM'V^d,), -g^deig. 

We will also denote by dvol = ^/gdx the natural Riemannian density, and (u, v) — Jj^j uvdvol 
the associated scalar product. 

In all the following, we will consider only the case where the waves are damped, wich 
corresponds to take a > with a non identically 0. We can reformulate the above problem 
into an equivalent one by considering the unbounded operator 

V -2ia ; 

with domain D{B) = H'^ x , and the following evolution equation : 

(1.3) ((9t + iS)u = 0, u = (uo,Ui) e i?^ X 
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From the Hille-Yosida theorem, one can show that B generates a uniformly bounded, strongly 
continuous semigroup e~'*^ for t > 0, mapping any {uq,ui) G x to a solution 
{u{t,x),idtu{t,x)) of (1.3). Since B has compact resolvent, its spectrum SpecB consist 
in a discrete sequence of eigenvalues {r„}„gN . The eigenspaces En corresponding to the 
eigenvalues r„ are all finite dimensional, and the sum 0„-B„ is dense in x H°, see 
[GoKr]. If r G SpecS, there is v £ such that 

(1.4) u{t,x)=e-'*''v{x), 
and the function u then satisfies 

(1.5) P{t)u = , where P(r) = - A - - 2 iar. 

From (1.5), it can be shown that the spectrum is symmetric with respect to the imaginary 
axis, and satisfies 

< ImT„ < 

while I ReT„| ^ oo as n ^ oo . Furthermore, if Rer ^ 0, we have Imr G [— ||a||oo, 0], and 
the only real eigenvalue is r = 0, associated to the constant solutions of (1.1). 

The question of an asymptotic density of modes has been adressed by Markus and Mat- 
saev in [MaMa], where they proved the following Weyl-type law, also found later indepen- 
dently by Sjostrand in [Sjo] : 

Card{n : < Rer„ < A} = ( f dxd£, + ©(A"^"^) . 

\27ry ip-i([o,i]) 

Here P — gx{£,,S,)'^ is the principal symbol of — Ag and dxd^ denotes the Liouville measure 
on T*M coming from its symplectic structure. Under the asumption of ergodicity for the 
geodesic fiow with respect to the Liouville measure, Sjostrand also showed that most of the 
eigenvalues concentrate on a line in the high-frequency limit. More precisely, he proved that 
given any £ > 0, 

(1.6) Card{n : r„ e [A, A -h 1] i(M \ [-a -e,-a + e])} = o(A'^"^) . 

The real number a is the ergodic mean of a on the unit cotangent bundle S*M = {(a;,^) € 
T*M,g4^,^) = 1}. It is given by 

a = ^lim J ao ^*dt, well defined dxd£, — almost everywhere on S*M . 

Hence the eigenvalues close to the real axis, say with imaginary parts in [a, 0], > a > —a 
can be considered as "exceptional". The first result we will present in this paper show that 
a spectral gap of finite width can actually exist below the real axis under some dynamical 
hypotheses, see Theorem 1 below. 

The second object studied in this work is the energy of the waves. From now on, we 
call H = X the space of Cauchy data. Let m be a solution of (1.1) with initial data 
u> GH. The energy of u is defined by 

E{u,t) = ^-mu\\h + \\yu\\h). 

As a well known fact, E is decreasing in time, and E{u,t) 0. It is then natural to ask 
if a particular rate of decay of the energy can be identified. Let s > be a positive number, 
and define the Hibert space 

= xH^ cn. 
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Generalizing slightly a definition of Lebeau, we introduce tfie best exponential rate of decay 
with respect to || • \\u^ as 

(1.7) p{s) = sup{/? G M+ : 3C> such that G W. E{uJ) < C c'^* \Mh'} 

where the solutions u of (1.1) have been identified with the Cauchy data u) G THf . It is 
shown in [Leb] that 

p(0) = 2min(G,C(oo)), 
where G = inf{— Imr; r G SpecB \ {0}} is the spectral gap, and 



Here : T*M — » T*M is the geodesic flow, and tt : T*M — » M is the canonical projection 
along the flbers. It follows that the presence of a spectral gap below the real axis is of 
significative importance in the study of the energy decay. However, an explicit example 
is given in [Leb], where G > while C(oo) = 0, and then p(0) = . This particular 
situation is due to the failure of the geometrical control, namely, the existence of orbits of 
the geodesic flow not meeting suppa (which implies C(oo) = 0). Hence, the spectrum of ,8 
may not always control the energy decay, and some dynamical assumptions on the geodesic 
flow are required if we want to solve positively the stabilization problem. In the case where 
geometric control holds [RaTa], it has been shown in various settings that p{i)) > 0, see for 
instance [BLR, Leb, Hit]. In [Chr], a particular situation is analyzed where the geometric 
control does not hold near a closed hyperbolic orbit of the geodesic flow: in this case, there 
is a sub-exponential decay of the energy with respect to || • ||-^£ for some e > 0. 

Dynamical assumptions. In this paper, we first assume (Af , g) has strictly negative 
sectional curvatures. This implies that the geodesic flow has the Anosov property on every 
energy layer, see Section 2.1 below. Without loss of generality, we suppose that the injec- 
tivity radius satisfies r > 2. Then, wc drop the geometric control assumption, and replace it 
with a dynamical hypothese involving the topological pressure of the geodesic flow on 5**-^, 
which we deflne now. For ervery s > and T > 0, a set 5 C 5*M is (e, T)— separated if 
p,9 <E S implies that (i((l)*/9, $*6') > e for some f G [0,T], where d is the distance induced 
from the adapted metric on T*M . For / continuous on S*M, set 



The topological pressure Pr(/) of the function / with respect to the geodesic flow is deflned 



The pressure Pr(/) contains useful information on the Birkhoff averages of / and the com- 
plexity of the geodesic flow, see for instance [Wal] for a general introduction and further 
properties. The particular function we will deal with is given by 

1 

(1.8) a^- ■.p&S*M^a^{p) = - 7r*a o $«(p) rfs - log J"(p) g M 

JO 2 

where J^{p) is the unstable Jacobian at p for time 1, see Section 2.1. In this paper, we will 
always assume that 



C(oo) = lim 





by 



Pr(/) = lim^ lim sup ^ log Z(/, T, e). 



(1.9) 



Pr(a") < 0. 



4 



EMMANUEL SCHENCK 



Main results. Under the condition Pr(a") < 0, we will see that a spectral gap of finite 

width exists below the real axis. As a consequence, there is an exponential decay of the 
energy of the waves with respect to || • Wn" for any k > d/2, and if G < | Pr(a")|, we have 
P{k) = 2G. We begin by stating the result concerning the spectral gap. 

Theorem 1. (Spectral gap) Suppose that the topological pressure of a" with respect to 
the geodesic flow on S*M satisfies Pr(a") < 0, and let e > be such that 

Pr(a") + £ < 0. 

Then, there exisits eo(£) > such that for any r S SpecB with | Rer| > eo(£), we have 

Imr < Pr(a") +£. 

The presence of a spectral gap of finite width below the real axis is not obvious a priori 
if geometric control does not hold, since there may be a possibility for |Imr„| to become 
arbitrary small as n ^ oo : see for instance [Hit], Theorem 1.3. However, this accumulation 
on the real axis can not occur faster than a fixed exponential rate, as it was shown in [Leb] 
that 

3C > such that Vr e SpccS, Imr < -l^c-C\Ker\ _ 

Let us mention a result comparable to Theorem 1 in the framework of chaotic scattering 
obtained recently by Nonnenmacher and Zworski [NoZw], in the semiclassical setting. For a 
large class of Hamiltonians, including P{K) = —hA + V on R'' with V compactly supported, 
they were able to show a resonance-free region near the energy E: 

3(5, 7 > such that Rcs(P(fi)) n{[E-6,E + d]- i[0, -fh]) = for < /i < hs,^. 

This holds provided that the hamiltonian fiow on the trapped set Ke at energy E is 
hyperbolic, and that the pressure of the unstable Jacobian with respect to the geodesic fiow 
on Ke is strictly negative. We will adapt several techniques of [NoZw] to prove Theorem 1, 
some of them coming back to [Anal, AnNo]. 

In a recent paper, Anantharaman [Ana2] studied the spectral deviations of Spec B with 
respect to the line of accumulation Imz = —a appearing in (1.6). In the case of constant 
negative curvature, she obtained an upper bound for the number of modes with imaginary 
parts above —a , and showed that for a G [— a, 0[, there exists a function H{a) such that 
(1.10) 

w ^ r, w ^ n 1- logCard{r„ : Rer„ e [A-c,A + c], ImT„ > a + £} 

Vc > 0,V£ > 0, limsup ; — r < H[a). 

A^oo log A 

H{a) is a dynamical quantity defined by 

H{a) ~ sup{/iifs(/i), /K G -Ml, J ad/i = —a} 

where Aii denotes the set of invariant measures on S*M, and hxs stands for the 
Kolmogorov-Sinai entropy of /x. As a consequence of Theorem 1, the result of Anantharaman 
is not always optimal : H{a) 7^ for a G [— a,0[, but if Pr(a") < 0, there is no spectrum 
in a strip of finite width below the real axis, i.e. the limsup in (1.10) vanishes for some 
a = a{a) ^ 0. 

The operator B being non-selfadjoint, its eigenfunctions may fail to form a Riesz basis 
of Ti.. However, if a solution u of (1.3) has initial data sufficiently regular, it is still pos- 
sible to expand it on eigenfunctions which eigenmodes are close to the real axis, up to an 
exponentially small error in time : 
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Theorem 2. (Eigenvalues expansion) Let e > such that Pr(a") + e < 0, and k> ^. 

There exists eo{e) > 0, n = n{e) € N and a (finite) sequence tq, ■ ■ ■ ,t„_i of eigenvalues of 
B with 

Tj e [-eo(e), eo(e)] + i[Pr(a") + £, 0], j e |0, n-l}, 
such that for any solution u{t,x) of (1.3) with initial data (jJ € W*, we have 

n-1 

u{t,x) = ^e"'*"^' Uj{t,x) +rn{t,x) , t>0. 

The functions Uj , rn satisfy 

\\uj{t, ■)\\n < Ct^' \Mh and ||r„(t, .)\\n < e'^''<"''^+'^ ||u;||„., 

where mj denotes the multiplicity of Tj, the constants C > depends only on M and a, 
while Ce > depending on M,a and e. 

A similar eigenvalues expansion can be found in [Hit], where no particular assumption 
on the curvature of M is made, however the geometric control must hold. Our last result 
deals with an exponential decay of the energy, which will be derived a consequence of the 
preceding theorem : 

Theorem 3. (Exponential energy decay) Let e > 0, (rj)o<j<„(e)_i, K and u as in 
Theorem 2. Set by convention tq = 0. The energy E(u,t) satisfies 

(\ 2 
n-1 

where rrij denotes the multiplicity of Tj . The constants C > depends only on M and a, 
while Ce > depending on M. a. and e. In particular, p{k) — 2min(G, | Pr(a") + e|) > 0. 

Remark. In our setting, it may happen that geometric control does not hold, while Pr(a") < 
0. In this particular situation, it follows from [BLR] that we can not have an exponential 
energy decay uniformly for all Cauchy data in H, where by uniform we mean that the 
constant C appearing in (1.7) does not depend on u. However, if for k > | we look at 
P{k) instead of p(0), our results show that we still have uniform exponential decay, namely 
p{k) > while p{0) = 0. 

1.1. Semiclassical reduction. The main stop yielding to Theorem (1) is more easily 
achieved when working in a semiclassical setting. From the eigenvalue equation (1.5), we 
are lead to study the equation 

P{t)u = 

where Imr = 0(1). To obtain a spectral gap below the real axis, we are lead to study 
eigenvalues with arbitrary large real parts since SpecS is discrete. For this purpose, we 
introduce a semiclassical parameter h e]0, 1], and write the eigenvalues as 

r=l+0(l). 

If we let /i go to 0, the eigenvalues r we are interested in then satisfy tH iizl£+ i Putting 
r = ^ and z = we rewrite the stationary equation 

^-^^ - 0-1%^^ M = 0, q;,{x) = V2za{x). 
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Equivalently, we write 

(1.11) {V{z,h)- z)u = 

where V{z,h) = — — i%z. The parameter z plays the role of a complex eigenvalue of 
the non-sclfadjoint quantum Hamiltonian V. It is close to the "energy" E = 1/2, while Imz 
is of order h and represents the "decay rate" of the mode. In order to recall these properties, 
we will often write 

(1.12) z^^ + hC, Ce C and Id =0(1). 

In most of the following, we will deal with the semiclassical analysis of the non-selfadjoint 
Schrodinger operator 'P{z, h) and the associated Schrodinger equation 

(1.13) ihdt'b = P{z,h)^ with H-J-H^s = 1. 

The basic facts and notations we will use from semiclassical analysis are recalled in Appendix 
A. The operator V has a principal symbol equal to p{x,(,) = ^gx{^,£,) , and a subprincipal 
symbol given by —iqz- Note that the classical Hamiltonian p{x,^) generates the geodesic 
flow on the energy surface p^^(^) = S*M. The properties of the geodesic flow on S*M which 
will be useful to us are summarized in the next section, where is also given an alternative 
definition of the topological pressure more adapted to our purposes. We will denote the 
quantum propagator by 

so that if * e L^(M) satisfies (1.13), we have *(t) = ZY**(0). Using standard methods 
of semiclassical analysis, one can show that is a Fourier integral operator (see [EvZw], 
chapter 10) associated with the symplectic diffeomorphism given by the geodesic fiow 
Since we assumed that a > 0, it is true that ||^/*||i,2_>j;,2 < 1, Vt > 0. 
Denote 

s 1 = = i + Oih) e C, 3* e l^{m), {v{z, h) - z)-<^ = 0}. 

2 2 

If ^; e El and ^ is such that (1.11) holds, the semiclassical wave front set of ^ satisfies 

WFr(*) c s*m. 

This comes from the fact that '5 is an eigenfunction associated with the eigenvalue ^ of a 
pseudodifferential operator with principal symbol p{x,S) = \gx{^^i)- Using these semiclas- 
sical settings, we will show the following key result : 

Theorem 4. Let ^; € S i , and e > be such that Pr(a") + £ < 0. There exists ho = ho{e) 
such that 

h<fki ^ — £<Pr(a") + £. 

h 

From (1.12), we also notice that the above equation impHes Imr < Pr(a") + e + 0{h) 
since r = h~^^/2z, and then Imr = Im^ + 0{h). It follows by rescaling that Theorem 4 is 
equivalent to Theorem 1. 

2. Quantum dynamics and spectral gap 
2.1. Hyperbolic flow and topological pressure. We call 

= e*^" : T*M T*M 
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the geodesic flow, where Hp is the Hamilton vector fleld of p. In local coordinates, 

rr defy^ dp dp 

Hp = 1^9il^^^-9^9,, = {p,-} 
1=1 

where the last equality refers to the Poisson bracket with respect to the canonical symplectic 
form uj = X]f=i '^^i ^ '^'^'i- Since M has strictly negative curvature, the flow generated by 
Hp on constant energy layers £ = p~^{E) C T*M, E > has the Anosov property: for any 
p € £, the tangent space Tp£ splits into flow, stable and unstable subspaces 

Tp£ = RHp®E%p)®E''{p). 

The spaces E^{p) and E"{p) are d — 1 dimensional, and are stable under the flow map: 

VteK, d^l{E%p)) = E%^\p)), d^l{E^{p))=E^{^\p)). 

Moreover, there exist C, A > such that 

i) \\d^l{v)\\ < C e-^* \\v\\, for all V e E^{p), t>0 

(2.1) ii) \\d^-*{v)\\ < Ce"^* \\v\\, for ah v e t > 0. 

One can show that there exist a metric on T*M call the adapted metric, for which one can 
takes C = 1 in the preceding equations. At each point p, the spaces E'^{p) are tangent to 
the unstable manifold W^{p), the set of points p" G f such that d($*(/9"), $*(p)) *^~°°) 
where d is the distance induced from the adapted metric. Similarly, E^{p) is tangent to the 
stable manifold the set of points such that (p*), $*((0)) *^+°°> q. 

The adapted metric induces a the volum form flp on any d dimensional subspace of 
T{T*M). Using fip, we now define the unstable Jacobian at p for time t. Let us define the 
weak-stable and weak-unstable subspaces at p by 

E''°{p) = E'{p) e Mi/p , = E''{p) e ^Hp. 

We set 

J'ip) = detrf$ * ij^.o($*(p)) = — — , J" p = Ji" p , 

where (v\, . . . , Vd) can be any basis of i?"'°(p). While we do not necessarily have J^{p) < 1, 
it is true that J"(p) decays exponentially as t ^ +cx). 

The definition of the topological pressure of the geodesic fiow given in the introduction, 
although quite straighforward to state, is not reafly suitable for our purposes. The alterna- 
tive definition of the pressure we will work with is based on refined covers of S*M, and can 
be stated as follows. For 5 > 0, let £^ = P^li ~ '^5 ^ + be a thin neighbourhood of the 
constant energy surface p~^{^), and V ~ {K>}ae/ an open cover of In what follows, we 
shall always choose 5 < 1/2. For T £ N*, we define the refined cover V'^'^\ made of the sets 

T-1 

F^=f|$-^(HJ, P = boh...bT-lGl^. 

It will be useful to coarse-grain any continuous function / on £^ with respect to y(^) by 
setting 

T-1 

(/)t,/3 = sup V/o$^(p). 
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One then define 

Zt{V, /) = mf J ^ exp((/)T,0) : Bt C , C \J vA . 

^ (l3eBT i3eBT J 

The topological pressure of / with respect to the geodesic flow on £^ is deflned by : 

Pr^(/)= lim lim llogZrCF,/). 

diam V^O 1 — ^oo 1 

The pressure on the unit tangent bundle S*M is simply obtained by continuity, taking the 
limit Pr(/) = lim^^o Pr*^!/)- To make the above limits easier to work with, we now take 
/ = a" and flx e > such that Pr(a")+e < 0. Then, we choose the width of the energy layer 
5 G]0, 1[ sufRciently small such that | Pr(a") - Pr'^(a")| < e/2. Given a cover V = {V„}ae^ 
(of arbitrary small diameter) , there exist a time to depending on the cover V such that 



ilogZt„(V,a-)-Pr^(a-) 

10 



e 

< - 
- 2 



Hence there is a subset of fo— strings Bt^ C such that {VajaeBtg is an open cover of £^ 
and satisfies 

(2.2) Y. exp((a")to,/3) < exp (io(Pr^(a") + |)) < exp (io(Pr(a«) + e)) . 

For convenience, we denote by \y^i3}i3eBt„ = {l^/3}/3eBto the sub-cover of V*-*"-* such that 

(2.2) holds. Note that in this case, the diameter of V, and then W depends on e. 

2.2. Discrete time evolution. Let {(/S/sj/JeBtu be a partition of unity adapted to W, so 

def 

that its Weyl quantization ip^ = 11^ (see Appendix A) satisfy 

WFR(n^) C £:^ = n^, ^ n^j = l microlocally near E^/"^ . 



We will also consider a partition of unity {<f>a}(xeA adapted to the cover V, and its Weyl 
quantization ft (p^ . In what follows, we will be interested in the propagator ^^*o+i ^ 
and 

N = Tlogh-\ T > 0. 

It is important to note that T can be arbitrary large, but is fixed with respect to h. The 
propagator U^*^" is decomposed by inserting ^Z/aeij^ each time step of length to- 

Setting first Uff = W*°n^, we have (microlocally near £^^^) the equality = and 
then 

(2.3) U^*" = J2 Ui3^...U0,, nearf'^/^ 

/3i,/32,.-,/3iveBto 

2.3. Proof of Theorem 4. We begin by choosing x & C§°{T*M) such that suppx € £^ 
and X = 1 on S''^^, and considering Op^(x). Applying the Cauchy-Schwartz inequality, we 
get immediately 

(2.4) 112^^*0+1 Op,(x) II < E m^3.■■■i^0^^'OMx)\\+OL^^L^{n'^). 

/3i,/32,...,/3NeB~, 
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Unless otherwise stated, the norms || • || always refer to || • ||i,2_>j,2 or || • \\l2, according to 
the context. The proof of Theorem (4) relies on the following intermediate result, proven 
much later in Section 4. 

Proposition 5. (Hyperbolic dispersion estimate) Let e > 0, and 6,V,to be as in 
Section 2.1. For N = Tto\ogh~^ , T > 0, take a sequence Pi, . . . Pn and Wi, . . . , the 

associated open sets of the refined cover W. Finally , let Op^(x) be as above. There exists 
a constant C > and ho{s) g]0, 1[ such that 

N 

h<ho ^ Wu'^Up, . ..u'°U0,u' Opft(x)|| < Ch-''/^ JJ e<»">*o.^, 

where {0^)10,0 = sup^gyy^ ^*j=o ° ^"'(p)- '^^^ constant C only depends on the manifold 
M. 

Wc also state the following crucial consequence : 

Corollary 6. Take £ > such that Pr(a") + e < 0. There exists C > and ho{e) G]0, 1[ 
such that 

h<ho^ + l Ops(x)|| < Ch-i eNto(Pr{an+e) 

The constant C only depends on M. 

Proof. Given e > 0, we choose 6,V,to,W as in the preceding proposition. Using (2.4), we 
then have 

ll^^'^*''+'Op,(x)|| < Ch-i J2 [lle^''"^^o.^^ +0{hn 

/?i.../3jveB,^ 

< Ch-i I c<"">'o.0| +0(fi°°). 

To get the second line, notice that the number of terms in the sum is of order (Card;8to)^ = 
^-TiogCardBtQ From our choice of e and S, we can use (2.2), and for h small enough, rewrite 
this equation as 

||ZY^*''+iOpft(x)|| < c;i-^e^*«(P^'('^")+^/2) 

< Ch-i gJVto(Pr(a")+£) 

where C > only depends on the manifold M. □ 



Let us show how this result imphes Theorem 4. We assume that ^ satisfies (1.11), and 

~ . Notice also that w( 

Opft(x)* = * + o{h- 



therefore = e* "^n' . Notice also that we have 



since WFr(*) C S*M, and then 

ll^Jvto+i Opft(x)*|| = \\U^*°+^^\\ + 0{h°°) = e^^*°+^^^'^'+0{h°°) . 
It follows from the corollary that 

g W < Ch-i e^*o(Pr(,")+e) ^ 
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where m can be arbitrary large. Taking the logarithm, this yields to 

But given e > 0, we can take A'' = Tlogfi"^ with T arbitrary. Hence there is ho{e) e]0, 1[ 
and T sufSciently large, such that 

Im z 

h < He) ^ — < Pr(g") + 2e. 
Since the parameter s can be chosen as small as wished, this proves Theorem 4. 

3. Eigenvalues expansion and energy decay 

3.1. Resolvent estimates. To show the exponential decay of the energy, we follow a stan- 
dard route from resolvent estimates in a strip around the real axis. Let us denote 

Q{z, h) = — -l^-z-i hV2za{x) = Viz, h) - z. 

The following proposition establish a resolvent estimate in a strip of width | Pr(a") + e| 
below the real axis in the semiclassical limit. This is the main step toward Theorems 2 and 
3, see also [NoZw2] for comparable resolvent estimates in the chaotic scattering situation : 

Proposition 7. Let e > 0. Choose 7 < such that 

Pr(a") + e < 7 < 0, 

and z = \ + hC,, with \C,\ = 0(1) satisfying 

7 < Im C < . 

There exists ho{s) > 0, Cg > depending on M, a and e such that 

h<ho{e) ^ ||g(z,n)-i|U2_i2 <Cefi-^+"'''"^log?i-^ 

where Co = 2|p,(g.)+,| - 

Proof. Given e > 0, we fix ?iq{£) so that Corollary 6 holds. Finally, we define 

Pr(a")+ =Pr(a")+e. 

In order to bound (5(z,/i)^^, wc proceed in two steps, by finding two operators which 
approximate Q^^: one on the energy surface £^ , the other outside £^ . Let x be as in 
Section 2, and choose also Let x G C^{T*M) with suppx suppx, such that we also have 
X = 1 near S*M. We first look for an operator ^0 — ^oi^, K) such that 

OAo = (l-Op;,(x)) + OL^^L^(ft°°). 

For this, consider (5(2;, ^) -|-iOp;j(x) =' Qo{z,h). Because of the property of x, the operator 
Qo is elliptic. Hence, there is an operator Aq, uniformly bounded in i^(M), such that 

The operator Aq we are looking for is obtained by taking Aq = Aq[1 — Opf((x)). Indeed, 
Q{z, h)Ao{z, h) = I- Opft(x) - iOpft(x)io(l - Opft(x)) + O^.^l^K^) 

since x and 1 — X have disjoints supports by construction. 
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We now look for the solution on £^ . From Corollary 6, we have an exponential decay of 
the propagator U^*° if N becomes large. To use this information, we set 

Ar{z, h, Tn) = ^ W e^i' Op^{x) dt 

where Tn has to be adjusted. Hence, 

Q(z, ft)Ai(z, ft, n) = Op;,(x) - W e^^ |t=T, Op^(x) = Op^(x) + Ri • 
Since Im z/ft = Im ^ = 0(1), we have 

(3.1) ||i?il|=e-^-i-«||Z^^-Op,(x)l|. 
From Corollary 6, we know that for A'' = Tig logft"^ with ft < fto, 

Opr,(x)|| < Ch-i e^*«P'-(«'')+ . 
We observe that this bound is useful only if does not diverge as ft ^ 0, which is the case if 

rp ^ d def 

- 2to|Pr(a»)+| ~ °' 

Let us define 

(3-2) n = Totologh-' + 1 = ^pj^^ log 71-^ + 1 • 

If Tfi = T\ogh~^to + 1, with T > To chosen large enough, we find 

< Cft^*"'^fi-^-"°P''(""'^ = ©(ft") 

with m = m{Tfi) > 0, since 7 — Pr(a")"'" > 0. Consequently, there is Ti = Ti{s) > such 
that = Tiiologft"^ + 1 satisfies m(T^) = 0. This means that for Tft > T^, we have 

Q{z,n){Ao{z,h) + Ai{z,h,Th)) = l + Oi2^i2(l), 

in other words, Aq + A\ is "close" to the resolvent . Hence, we impose now T% > T^, 
and evaluate the norms of ^0 and Ai. By construction, ||j4o|| = 0{1). For Ai, we have to 
estimate an integral of the form 

lTn= ||iY*Op,(x)||rft. 
^0 

Let us split the integral according to T^, and use the decay ofU* Op^(x) for t>T^: 

nOO 

\It,\ < rOe-^^'"«+Cft-* / e-"'"«e(*-i)P'-(«'')+rft 

< TOe-^"i-«(l + C,h-i e(^-°-i)P--(«")^) 

< Cer°e-^^°i'"« . 

Using (3.2), this gives 

Pi(^,^,Tr)|1 <aft-i+^°i'"^logft-i 
where Cg > depends now on M, a and £ while 

^'■'^ = 2|Pr(l«)+r 

□ 
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We now translate these results obtained in the semiclassical settings in terms of r. Recall 
P(r) = - A - - 2 iar = ^Q{z, h). 

and set R{t) =^ P(t)^^. The operator R{t) is directly related to the resolvent (r — B)~^ : 
a straightforward computation shows that 

MW(-2ia-T) -R{t)\ 
~ \ R{T)i2iaT-T^) -R{t)t )■ 

Proposition 8. Let e > be such that Pr(a") + e < 0, and 7 < satisfying 

Pr(a") +e < 7 < 0. 

Let r G C \ SpecS be such that 7 < Imr < 0. Set also (r) = (1 + |Tp)2 . There exists a 
constant C > depending on M, a and e such that for any k > d/2, we have 

(i) II i?(T) |U2^^. <Ce(r)-i-^»^log(r) 

(^Z) ||i?(T)|U2^^2 <Ce(T)i-<=<'^log(T) 

(ill) \\R{t)\\h''^H^ <Ce 

(iv) \\TR{T)\\Hr^^HO <Ce. 

Proof, {i) fohows directly from rescahng the statements of the preceding proposition. For 
{ii), observe that 

\\RiT)u\\H'. < C{\\R{T)uh'2 + \\AR{T)uh2) ,00. 

But 

||Aii(T)M||z,2 < \\u\\l2 + + 2Tia\\\R{T)u\\L2 , 

so using {i), we get 

\\R{t)u\\h2 <C{{1 + \t^ + 2iaT\)\\R{T)u\\L2 + \\u\\l2) < Ce(r)i-'=°Tlog(T)||w|U2 . 
To arrive at (Hi), we start from the following classical consequence of the Holder inequality 

(3.4) \\R{t)u\\%,-. < \\R{T)uf-:>^R{T)u\\]+' , s>0. 

From {i) and {ii), we obtain 

\\R{T)u\\m-s < C,{T)-^^'°+'hog{T)\\u\\L2 . 

If we choose s > — 7C0, we get \\R{T)\\f[o^fji-s < C^- Hence, for any s' > we have 

\\R{T)\\H.'^H'' + ^-'<Ce- 

Taking s' = s shows {Hi), where we must have k > —jcq. In view of (3.3), and the fact 
that 7 > Pr(a^) + e, this condition is satisfied as soon as k > d/2. The last equation {iv) is 
derived as {Hi), by considering 

\\rR{r)u\\l.-. < \TnR{r)u\\],r\\R{r)u\\]+^ , s > 0, 

and choosing s so that ||T-R(r)||^o_>jifi-s < C'e- D 
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3.2. Eigenvalues expansion. We now prove Theorem 2. Let us fix e > so that Pr(a") + 
£ < 0. From Theorem 1 we know that 

Card (SpecB n (R + i[Pr(a") + e, 0])) =^ n{s) < oo. 

Hence there is eo{s) > such that SpecS n (M + i[Pr(a") + e]) C f2, where 

n = n{s) = heo, ea] + i[Pr(a'') + £, 0]. 

We then call {tq, . . . , r„(e)_i} = SpecB D fl, and set by convention tq = 0. We define as 
above Pr(a")+ = Pr(a") + e. Since we look at the eigenvalues r G fi, let us introduce the 
spectral projectors on the generalized eigenspace Ej for j S |0, n — Ij : 

where jj are small circles centered in Tj. We also denote by 

the spectral projection onto 0j=o^j- We call Eq the eigenspace corresponding to the 
eigenvalue tq = 0. It can be shown [Leb] that Eq is one dimensional over C and spanned by 
(1,0), so 

Hoo; = (c(a;),0) with c{uj) G C. 
Let now w = {lvo,loi) be in Ti'^. Near a pole tj of (r — B)~^, we have 

where Hj is an operator depending holomorphically on r in a neighbourhood of tj , and rrij is 
the multiplicity of Tj. Since H e jC{H, D{B°°)), we have the following integral representation 
of e~'*^na;, with absolute convergence in H: 

1 (•+oo+ia 

(3.5) e-'*^nu;=— / e-'*^(r - B)-^nwrfr , i > 0, a > 0. 

2 ITT ./-oo+i a 

The integrand in the right hand side has poles located at tj, j G |0,n — 1], so that 




The operators (t) appearing in the residues are polynomials in t, with degree at most 
rrij, taking their values in C{Ti., D{B°°)). It follows that for some C > depending only on 
M and a, 

lh,(t)||«<Ct"^||w||«. 
The remainder term appearing in Theorem 2 is now identified : 

r„(t)=e-'*'^(l-n)u;. 

To conclude the proof, we have therefore to evaluate ||r„||-K. To do so, we will use in a 
crucial way the resolvent bounds below the real axis that we have obtained in the preceding 
section. We consider the solution u{t,x) of (1.1) with initial data u = {uo,Ui) = (1 — II)w, 
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with ugW^, K> d/2. Let us define x e C°°(IR), < x < 1, such that x = for i < and 
X = 1 for t > 1. If we set u = x^, we have 

(3.6) {d^ - A + 2adt)v = gi 
where 

(3.7) gi = x"u + 2x'dcu + 2ax'u . 

Note also that supp^i C [0, 1] x M, and v{t) = for f < 0. Let us denote the inverse Fourier 
transform in time by 

:Ft^-r : u ^ u{t) = / c'*'" u{t)dt. 

Applying J-t^-r (in the distributional sense) to both sides of (3.6) yields to 

P{t)v{t,x) = gi{T,x) . 
We then remark that i?(r)gi(T, x) is the first component of 

i{T - B)-'j^t^-r {x'mu,idtU)). 

From the properties of 11, it is clear that the operator (r — B)~^{1 — TV) depends holomor- 
phically on r in the half-plane Imr > Pr(a")+. From {u,\diu) = c"'*®(l — n)a;, we then 
conclude that i(r — B)~^ J^t^-T {x! {i){u^\dtu)) depends also holomorphically on r in the 
half plane Imr > Pr(a")"'". Hence v{t,x) = R{T)g\{T,x) and an application of the Parseval 
formula yields to 

||e-*P''("")%(t,a;)|U2(M^,Hi) = ||i;(r + iPr(a'')+|U.(M,Hi) 

= \\R{t + iPr(a")+).gi(T + i Pr(a-)+, x)|U2(k,hi) 

< C,||.gi(T + iPr(a")+,a:)||i2(K,H'>) 

< 11,91 (t..x) II L2(R^,ff») . 

where we have used Proposition 8. The term appearing in the last line can in fact be 
controlled by the initial data. From (3.7), we have 

(3-8) ||5'l||l,2(M+;H'^) < C (||u||i2([0,l];/i-'') + || « || L2([0,l] iff-^ ) ) • 

A direct computation shows 

dtMl^. < c{\\u\\\i. + \\dtu\\]j. + \\vu\\]j.) . 

The Gronwall inequality for t S [0, 1] gives 

\\u{t,.)\\%. < c(^imo,-)||l,. + ^*(||ax«)||l^. + ||v^.)||l,.)rf^ 

< C-llu^ll?,., 

since the k— energy 

E-{t,u) = \murH^ + \\'^u\\i.) 

is also decreasing in t. Coming back to (3.8), we see that ||,9i||l2(r^.j:^^-) < CHwH-h- and 
then, 

||e-*(P'-('^")+^)z;(L:r)|U2(R^,^i) < a||a;||«.. 

This is the exponential decay we are looking for, but in the integrated form. It is now easy 
to see that 
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We have to check that the same property is valid for dtU. Using the same methods as above, 
we also have 

P{T):Ft-,-r{dtv) = ~Tgi{T), 
and then, Tt^^r{dtv) — ~TR{T)gi{T). It follows that 



Grouping the results, we see that 
and this concludes the proof of Theorem 2. 

3.3. Energy decay. We end this section with the proof of the Theorem 3, which gives the 
exponential energj' decay. This is an immediate consequence of the following lemma, that 
tells us that the energy can be controlled by the i? ^ norm of u, for t > 2: 

Lemma 9. There exists C > such that for any solution u of (1.1) and E{u,t) the 
associated energy functional, we have 

E(«,r)<C||^/||i.([^_2,T+i].^.), T>2. 

Proof. This is a standard result, we borrow the proof from [EvZw]. For T > 2, we choose 
X2 G C°°(R), < X2 < 1 such that X2(i) = 1 for t > T and X2{t) = if i < T - 1. Setting 
U2{t,x) = X2{t)u{t,x), we have 

{df - A + 2adt)u2 = 52 
for g2 = + 2x2^tW + 2ax2W. Note that g2 is compactly supported in t. Define now 

E2{u,t) = \ I {\dtU2? + |Vw2|^)rfvol 

and compute 



= ||i)(T + iPr(a")+|U2(K,ao) 

= \\tR{t + iPr(a")+).gi(r + iPr(a")+, x)|U2(r,^o) 

< C,||5i(T + iPr(a")+,x)|U2(K,^^») 

< Ce\\gi{t,x)\\L2i^^^^H-)- 



E'^iu, t) = {dfu2, dLU2) - {AU2, d,U2) 

= -2{adtU2, dtU2) + {g2,dtU2) 

< C [ \dtU2\{\dtU\ + \u\)dYO\ 

< c(E2{u,t)+ [ 



M 

i\dtu\^ + \u\^)dYO\ 



We remark that E2{u,T — 1) = and E2{u,T) = E{u,T), so the Gronwall inequality on 
the interval [T- 1,T] gives 

(3-9) E{u,T) < C (!|<9tM||i2([T-l,T]:L2) + ||w||i2([T-i,T];L2)) • 

To complete the proof, we need to bound the term ||5tw|||^2(['r_i 7'].j;,2). For this purpose, 
we choose xs S C°°(IR), < xs < 1 such that X3(0 = 1 for t e [T - 1,T] and X3(i) = if 
i < T - 2 and t > T + 1. From (1.1), we get 
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/ (X3^) — Au + 2adtu)dt 

Jt-2 
rT+1 

/ -xlidtu, dfu) - 2x3X3 (w, dtu) + 2xi(u, adtu) + xl{u, -Au)dt , 
Jt-2 



whence 



\\dtu\\L2([T-l,T];L^) < C'||M||L2([T-2,T+l];i/i)- 

Substituting this bound in (3.9) yields to the result. □ 

The Theorem 3 follows now from the preceding lemma. Let us denote by Uj {t, x) and and 
rn{t,x) the first component of Prj{t)u> and e~'*^(l — n)a; respectively. We learned above 
that 

n 

u{t, a;) = ^ c~ Uj{t, x) + rn{t, x) 

3=0 

with \\uj{t,-)\\Hi < Ct^^Mw, and ||r„(t,-)||j/i < e*(P'-('^")+2-) ||u;||„.. Suppose first 
that the projection of a; on Eq vanishes, i.e. Ilow = 0. Then, from the preceding lemma we 
clearly have 

n 

E{u,t)i < 5;^e"'"^^ C||«,(t,a;)||Hi +C,e*(P'-("")+2e) ||^||^^ _ 

This shows Theorem 3 when IIow = 0. But the general case follows easily: we can write 
u{t,x) = u{t.x) — IIow for which we have the expected exponential decay, and notice that 
E{u, t) ~ E{u, t) since now is constant. 



4. Hyperbolic dispersion estimate 

This last section is devoted to the proof of Proposition 5. Let e,5,V,W and Op^(x) be 
as in Section 2. We also set N = T loghr'^ , T > 0. 

4.1. Decomposition into elementary Lagrangian states. Recall that each set W/j = 
W'f,o...6t|j_i in the cover W has the property 

(4.1) $'=(W/3)cH,, fce[0,to-il 

for some sequence 6o, &i, • ■ • , ^to-i- We will say that a sequence sequence /3i . . . /3jv of sets 
W|3^. is adapted to the dynamics if the following condition is satisfied : 

Vfce[i,7V-i], <i>'=*''(>Vr,Jn>V/3,^, y^0. 

In this case, we can associate to the sequence {j3i} a sequence 7i, . . . ,7jvto of sets V^^ C V 
which are visited at the times 0, . . . , Nto—1 for some points of W/3i . We will only consider the 
sequences adapted to the dynamics. Indeed, it is clear from standard results on propagation 
of singularities that 

\\Up,...Up,\\=0{h^) 

if the sequence is not adapted (see Appendix A), and in this case. Proposition 5 is obviously 
true. 
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We now decompose further each evolution of length to in (2.3) by inserting additional 
quantum projoctors. To unify the notations, we define for j € |1, Ntoj the following projec- 
tors and the corresponding open sets in T*M : 

.42^ p ^i^P. ifj-l = fcto, /ceN ^ ^1^/3. iij-l = kto, keN 

[n^^. if j- 1^0 mod to ' \v^^ ifj-17^0 mod to- 

We will also denote by fj G Cq^{T*M) the function such that supp C V-^ and P-y = F™. 

Let us set up also a notation concerning the constants appearing in the various estimates 
we will deal with. Let LK ^ N be two parameters (independent of K), and 61,62,63 > 
some fixed numbers. For a constant C depending on M and derivatives of x, a, <&* (for t 
bounded) up to order ei^ + 62-?!' + 63, we will write C^^'-^\M,x)^ or simply C'^^\M,x) if 
only one parameter is involved. If the constant C depends also on the cutoff functions F^ 
and their derivatives, we will write 

C = C(^'^)(M,x,V). 

This is to recall us the dependence on the cutoff function x supported inside , and the 
refined cover V. We will sometimes use the notation C^^'^^(M, V) when no dependence on 
X is assumed. Note that V depends implicitely on e since its diameter was chosen such that 

(2.2) holds. 

Using (4.1), standard propagation estimates give 
and similar properties for U*"Ilpf^, k > 1. Finally, 

(4.3) U/s,,^ . . . U/s.U Opft(x) = WP^,,„ . . MP^^U OpM + Ol^^l^ {h^) . 

Take now '5 G L'^{M). In order to show Proposition 5, we will write Op^(x)'I' as a Hnear 
decomposition over some elementary Lagrangian states, and study the individual evolution 
of such elementary states by ZY^*""*"^. This type of method comes back to [Anal] and is the 
key tool to prove Proposition 5. The decomposition of Op;j(x)^' is obtained by expliciting 
the action of Op;j(x) in local coordinates (see Appendix A). When applying Op^(x) to \1/ 
using local charts labelled by we get 

[Op,(x)*](x) = X^^^|e'^^^x(^,r?)v'K^o)</'H^)*Mrfr;dzo 
where we have defined 

This is a Lagrangian state, which Lagrangian manifold is given by 

A° =^ T* MnS^ C T*M . 

Geometrically, A" corresponds to a small, connected piece taken out of the union of spheres 
{T*^M r\p~^{^ + y) , \v\ < 5}. If we project and evolve ^ according to the operator 
appearing in the right hand side of (4.3), we get : 
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£ z J M 

(4.4) < CY,^uv\\U'P^,,^...P^M5l.Mn 

where C > depends only on the manifold M. Hence we are lead by this superposition 
principle to study in detail states of the form U^P^^ . . MP-^-^^US^ .^^, for n G pjA^io] and 
t G [0, 1]. For simpHcity, because the local charts will not play any role in the following, we 
will omit them in the formula. 

4.2. Evolution of Lagrangian states and their Lagrangian manifolds. 

4.2.1. Ansatz for short times. In this section we investigate the first step of the sequence of 
projection-evolution given in (4.3): our goal is to describe the state L{ S^^zo 

with t e [0,1]. 

Since is a Fourier integral operator, we know that U^S-^^zg is a Lagrangian state, supported 
on the Lagrangian manifold 

AO(t) "li:^ $*(A0) , [0,1]. 

Because of our assumptions on the injectivity radius, the flow : A°(s) — > A°(t) for 
1 > t > s > 0, induces on M a bijection from 7rA°(s) to 7rA°(i). In other words, A°{t) 
projects diffeomorphically on M for t e]0, 1], i.e. kerd7r|A0(t) = : in this case, we will say 
that A*^(t) is projectible. This is the reason for introducing a first step of propagation during 
a time 1 : the Lagrangian manifold A"(0) is not projectible, but as soon as t s]0, 1], A^{t) 
projects diffeomorphically. Treating separately this evolution for times t € [0, 1] avoid some 
unnecessary technical complications. 

The remark above implies that the Lagrangian manifold A"(t), t e]0, 1] is generated by 
the graph of the differential of a smooth, well defined function 5*0 : 

A°{t) ={{x,da;So{t,x,zo)) : 1 > t > 0, x e 7r$*(A°)} . 

This means that for t e]0, 1], we have the Lagrangian Ansatz : 

v°{t,x,zo) U^S^^zoix) 

(4.5) = ^^[e'''''^'^hHUt,x,zo)+h''B'i,{t,x,zo)]. 

The functions b^{t,x,zo) are smooth, and x G 7rA°(t). Furthermore, given any multi index 
£, they satisfy 

(4.6) \\dibl{t,;Zo)\\<Ce,k 

where the constants Ce^k depends only on M (via the Hamiltonian flow oip), the damping a, 
the cutoff function x and their derivatives up to order 2k + £. However, note that Co.o only 
depends on M. The remainder satisfies H-Bxll — where the constant Ck also depends 
on M, a, X and is uniformly bounded with respect to x,zo. The base point zq will be fixed 
until section 4.5, so it will be ommited in the following to simplify the notations. 
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4.2.2. Further evolution. In the sequence of projection-evolution (4.3), we then have per- 
formed the first step, and obtained an Ansatz forU^S^, t e]0, 1] up to terms of order h^^''^^, 
for any K > 0. The main goal of the next paragraphs consist in finding an Ansatz for the 
full state 

(4.7) v"-{t,x) =U*P^JAP-y^_, . . .UP^JAS^ , t e [0, 1], n > 1 . 

The I3j are defined according to j — 1 mod as in the preceding section, but here n is 
arbitrary in the interval |l,iVfo]. Because the operator WP is a Fourier integral operator, 
v^{t,x), j > 1 is a Lagrangian state, with a Lagrangian manifold which will be denoted 
by A-'(f). This manifold consist in a small piece of (A°), because of the successive 
applications of the projectors between the evolution operator W. If j = 1, the Lagrangian 
manifold A^(0) is given by 

Ai(0) = AO(l)nV^„ 

and for t e [0, 1] we have A'^{t) = $*(Ai(0)). For j > 1, A^{t) can be obtained by a similar 
procedure: knowing A^^^(l), we take for A-'(t), t G [0, 1] the Lagrangian manifold 

A-''(0) =^AJ-i(l)nV^^ , and A^'(t) = $*(A^ (0)) . 

Of course, if the intersection A-'~^(l) fl is empty, the construction has to be stopped, 
since by standard propagation estimates, will be of order 0{h°°). But this situation will 
not happen since the sequence {Pi} is adapted to the dynamics. It follows that 

Vie[l,nl, A^(O)^0. 

One can show (sec [AnNo], Section 3.4.1 for an argument) that the Lagrangian manifolds 
A^{t) are projectible for all j > 1. This is mainly because M has no conjugate points. In 
particular, any A^{t) can be parametrized as a graph on M of a differential, which means 
that there is a generating function Sj{t, x) such that 

A^{t) = {x,d,S,(t,x)}. 

By extension, we will call a Lagrangian state projectible if its Lagrangian manifold is. 

Let us introduce now some notations that will be often used later. Suppose that x G 
nA^t), j > 1- Then, there is a unique y = y{x) G 7rA-'(0) such that 

no^*{y,dySj{0,y))=x. 

If we denote for t S [0, s] the (inverse) induced flow on M by 

^s*(s) • ^ ^ 7rA-'(s) 7r$~* {x, dxSj{s, x)) € 7rA-^(s - t), 

we have y{x) = <j)~^^, ^^■^{x) . If a; € 7rA-'(t), then by construction 

^-'-\x,d,Si{t,x)) e Ai-\Q) c A^-'=-i(i) , k e I0,i - 11 . 

By definition, we will write 

<PsXt)^^^ = -K^~*~''{x,dxSj{t,x)) and (i>s^{x) = '!r^~''{x,dxSj{l,x)) . 

To summarize, our sequence of projections and evolutions can be cast into the following 

way: 

(4.8) v°{l, .) v\0, •) v\l, •) ^ v^{0, •) v-{t, •) 

AO A0(1) Ai(0) A\l) ^ A"(0) A"(i) 
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On the top line are written the successive evolutions of the Lagrangian states, while the 
evolution of their respective Lagrangian manifolds is written below (the notation |v denotes 

a restriction to the set V C T*M). 

4.3. Evolut ion of a projectible Lagrangian state. Let and be as in (4.2). The 
next proposition contains an explicit description of the action of the Fourier integral oper- 
ators U^P on projectible Lagrangian states locahzed inside V^. 

Proposition 10. Let V-y and P-y = be as in Section 4-1- Let 'Wn{x) = w(x)e^'^(^' be a 
projectible Lagrangian state, supported on a projectible Lagrangian manifold 

A = {x,dxip{x)} C . 

Assume also that A{t) $*A is projectible for t € [0, 1]. We have the following asymptotic 
developement : 

K-\ 

(4.9) [U'P^wn\{x) = e^^(*'^) fi''Mt,x) + h''rK{t,x) 

fe=0 

where fl^it,-) is a generating function for A(t). The amplitudes Wk can be computed from 
the geodesic flow (via the function ^p-y), the damping q and the function Fj. Moreover, the 
following bounds hold : 

||w^fc|lc' < Ctk\\w\\ce+2k 
\\rK\\c' < C£,x||w|lc"!+2if+d 

where the constants depend on ip~i,, a, F-^, and their derivatives up to order £ + 2K + d, namely 
Ct^k = C^^''^''(M, V). An explicit expression for Wk will be given in the proof. 

Proof. The steps we will encounter below are very standard in the non-damping case, i.e. 
q = 0. If the diameter of the partition V of £^ is chosen small enough, we can assume 
without loss of generality the existence of a function ip../ G C°°([0, 1] x M'^ x R'^) which 
generates the canonical transformation given by the geodesic flow on for times t G [0, 1], 
in other words : 

(4.10) V(?/,77)eV^, ^*{y,r]) = {x,^)<^^ = da,ip^{t,x,r]) and y = dnipy{t,x,r]). 
Furthermore, ipj satisfies detd^ .^(pj ^ 0, and solves the following Hamilton-Jacobi equation 

(dtip^ +p{x,da:ip^) = 
\(Pj{0,x,r]) = {r],x) . 

We first look for an oscillatory integral representation: 

{4.np%wn{x) = jj e^('^-(*'-''')-fe'''>+^to)) ^ h'^alit, x, y, v)w{y)dydr, 

= bn{t,x) + n''fK{t,x), \\fK\\=0{l), 

with {y,ri) G . For simplicity, we will omit the dependence on 7 in the formulae. We 
have to determine the amplitudes a^. For this, we want bn to solve 
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up to order . Direct computations using (1.2) show that the functions ip and must 
satisfy the following equations : 

{dtif + p{x, dxip) = (Hamilton-Jacobi equation) 

dtao + qao + X[ao] + ^ao divg X = (0-th transport equation) 

dtcik + qak + X[ak] + divg X = ^Aga^-i {k — th transport equation) 
with initial conditions 

^{0,x,r]) = {x,r]) 
ao{0,x,y,7j) = F(^,J7) 
afc(0, X, y,ri) = for fc > 1. 

The variables y and rj are fixed in these equations, so they will play the role of parameters 
for the moment and will sometimes be skipped in the formulae. X is a vector field on M 
depending on t, and diVg X its Riemannian divergence. In local coordinates, 

X = g^'{x)dxjip{t,x) = d^^p{x,dx^p{t,x))dxi and divgX = ^9j(v^X'). 

v9 

The Hamilton-Jacobi equation is satisfied by construction. To deal with the transport 
equations, we notice that X corresponds to the projection on M of the Hamiltonian vector 
field Hp at {x,dx<f{t,x,rij) G T*M. Let us call first 

At.rj = {{x,dxipit,x,ri)),x e 7r$*A}, t] fixed. 

This Lagrangian manifold is the image of the Lagrangian manifold Aq,,, = {{y,il) '■ y & t^^} 
by the geodesic flow The flow «;* on M generated by X can be now identified with the 
geodesic flow restricted to Ag^^: 

K* : ttAs^t, B Xh^ '!r^*{x,dxip{t,x,r))) e TrAt+s,n- 

The inverse flow {kI)~^ will be denoted by k,7+c Let us extend now the flow k* of X on M 
to the flow /C* generated by the vector fleld X = dt + X on R x M : 

I'm X M -> M X M 
■ \ {s,x) 1-^ {s + t,Ki{x)). 

We then identify the functions Ofc with Riemannian half-densities on Mx M - see [Dui, EvZw]: 



ak{t,x) = ak{t,x)^ydtdvo\{x) = ak{t,x)^Jg{x)^ \dtdx\^ G C°°(M x M,ni). 
Since we have 

Cx{aWdtdvo\) = {X[ak] + -ak divg X)Vdtdvo\ , 
the 0-th transport equation takes the simple form of an ordinary differential equation: 



This is the same as 



which is solved by 



Cx{ao^dtdvoV) + qaoV dtdvol = 0. 



"T" (/C* ) * ao V dtdvol = —qaoV dtdvol, 
dt 



aoVdtdvol = e-/o''°'=° ''^%IC-*)* aoV dtdvol . 
We now have to make explicit the coordinates dependence, which yields to 

ao{t, x)^/§(^^\dxdt\ 5 = e" •''o 9°«r*(^)rf^ ^^(o, Kt\x))^ g{Kt\x))\ det d^K^*]^ \dxdt\ 5 . 
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Consequently, 

V9{x) 

Since 

: X 7r$~*(a;, a^yjCt, x, r/)) = dj^Lpit, x, rj), 

it is clear that | det = I det a;, ?7)|. For convenience, we introduce the fol- 

lowing operator transporting functions / on M with support inside TrAg,,, to functions 
on TrAt+s,rj while damping them along the trajectory : 

This operator plays a crucial role, since we have 

(4.13) oo(t,-)=To*(ao(0,-))=^oF, 

from which we see that ao{t,-) is supported inside ttAj^^. By the Duhamel formula, the 
higher order terms can now be computed, they are given by 



ak{t,-) = 1^ T/-« (^iAgOfc_i(s)^ ds. 



The ansatz bh{t,x) constructed so far satisfies the approximate equation 

at 

The difference with the actual solution U^P is bounded by 

h^t\\A„aK-i\\ < Cth^ , 



= (i hAg - q)bn -^^"^ II e^^^*'^'"'^) w{y)AgaK-i{t, x, y, v) dydr, . 



where C = C(2^)(M, V), so (4.11) is satisfied. 

As noticed above, for time t > 0, the state U^Pwh. is a Lagrangian state, supported on 
the Lagrangian manifold A{t) = $*A. By hypothese, A(t) is projcctiblc, so wc expect an 
asymptotic expansion for bfi(t,x), exactly as in (4.5). To this end, we now proceed to the 
stationary phase developement of the oscillatory integral in (4.11). We set 

^'^""^ " J2^ II ^^^^^'''''^^~^'''^^^^^'^^Mt,x,y,v)w{y)dydv. 
The stationary points of the phase are given by 

f^'(y) = V 

\dn(p{t,x,ri) = y, 

for which there exists a solution (jjcilc) S A(0) in view of (4.10). Moreover, this solution 
is unique since A{t) is projectible: yc = yc(x) G 7rA(0) is the unique point in 7rA(0) such 
that X = n^*{yc,i>'{yc)) , and then rjc = V''(2/c) is the unique vector allowing the point yc to 
reach x in time t. The generating function for A{t) we are looking for is then given by 

ip{t, x) = S{t, X, ycXx), Vcix)). 
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Applying now the stationary phase theorem for each Ik (see for instance [Hor], Theorem 
7.7.6, or [NoZw], Lomma 4.1 for a similar computation), summing up the results and ordering 
the different terms according to their associated power of h, we see that (4.14) holds with 

|det(l - d^^ip{t,x,i^c) o (yc))h 

and 

k 

(4.14) Wk{t,x) = ^A2i{x,D^^r,){ak-i{t,x,y,r])w{y))\^y^ri)={yc,vc) ■ 

denotes a differential operator of order 2i, with coefficients depending smoothly on (p, 
V' and their derivatives up to order 2i + 2. This yields to the following bounds : 

||w^fe||c* < C'^,fe||w^||c^+2'= 

where Ce^k = C^^''^^(M, V). The remainder terms rKit,x) is the sum of the remainders 

coming from the stationary phase developement of Ik up to order K — k. Each remainder 
of order K — k has a C^' norm bounded by C^^^-fc/i^^'^ll w||(^€+2(/f-fc)+d so we see that 



< CiM\\w\\ci+2K+^, C = C^^'^\m,V). 



The principal symbol wq can also be interpreted more geometrically. As in Section 4.2, 
denote by 4>^^(^t) following map 



-t 



7rA(t) ^ 7rA(0) 



' \x^TT^-\x,d^ij{t,x)). 

Let us write the differential of $* : (t/,ry) <—^ (a;,^) as d^*{Sy,6r]) = {Sx,S^). Using (4.10), 
we have 



and then, since d^^^Lp is invertible 



If we restrict to A(0), we have Sij = xp"{y)Sy, which means that for x € 7rA(f), 

'^'t'i,it)(x) = ^^r,'^(*>2;,77c)(l - d^^(p{t,x,r]c)^p"{yc))~^. 
It follows from (4.14) that 



wo{t,x) = e^''(*)«;(ye)F(2/e,r?c) e" ^ 1 det #^[,)(x)| ^ V 

= ei/'W w{yc)f{yc, Vc) e" ^° «(^^«^(-))<^^ | Jac(#;*,,(a;)) 

where Jac(/) denotes the Jacobian of / : M — > M measured with respect to the Riemannian 
volume. □ 
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4.4. Ansatz for n > 1. In this paragraph, we construct by induction on n a Lagrangian 
state b"'{t,x) supported on A"(<;), in order to approximate v'^{t,x) up to order h^~'^^^. 

Proposition 11. There exists a sequence of functions 

{bl{t,x),Sn{t,x) : n>l, k < K, x € M, te[0,l]} 

such that S„{t,x) is a generating function for A"(t), and 

(4.15) v"{t,x) = ^^e'^ ^ h''bUt,x)+h''-^R'k{t,x) 

where BJ^ satisfies 

Cn K-1 \ 
EEii^r'(i,-)iic^(--m+c' 
i=2 k=0 / 

where C = C^^^Ad, x), Ck = C(^)(Af, x, V) and C > Q is fixed. 

Proof. The construction of the amplitudes 6^' for all A: > is done by induction on n, 
following step by step the sequence (4.8). In Section 4.2 we obtained U^S^ as a projectible 
Lagrangian state: 



k=0 



and we know that 6°(1,-) satisfies the hypotheses of Proposition 10, which will be used to 
describe U*P.y^v^ {!,■). More generally, suppose that the preceding step has lead for some 
n > 1 to 

(27r/l)2 ^ 
= -\-Tb''-Ht,x) + h^-'/^Rl-\t,x) 

where b"~^{t, •) is a Lagrangian state, supported on the Lagrangian manifold A"~^(t), and 
is some remainder in L^(M). We now apply Proposition 10 to each Lagrangian state 
Q-^Srv-i(i,x) appearing in the definition of 6"^^. Because of the term ft*', if we 

want an Ansatz as in (4.15), it is enough to describe U*'P^^v''^^^ {1, •) up to order K — k, 
which gives a remainder of order CK-kff'^''''Wk^'^{'\-,-)\\cnK-h)+d. Grouping the terms 
corresponding to the same power of h when applying Proposition 10 to each {v^~^)Q<k<K-i 
yields to 

K-l 

[W*P^„6"-i](x) =e^^"(*'^) E hHl{t,x)+h^B'}^{t,x) b''{t,x)+h'^B'^{t,x), 

fe=0 

where Sn{t,x) is a generating function of the Lagrangian manifold 

A"(i) = $*(A"-i(i)nv^j. 
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The coefficients 6^ are given by 

k k — i 

(4.17) 6£(t,a;) =^^^2z(a;j",_,(t,x,y,7?)6r'(l,2/))l(y,r,)=(j,.,r,.) 



i=0 1=0 

where = (j)g*^^^{x) , rjc = dyS„-i{l,yc). In particular, W${t,x) = D„(t, a;)6o~^(l, i/c), with 

(4.18) I?„(i,x) - c- J"" (^»''-^- |jac(d,^5^(,)(x))|%*'^"W F^„(ye,r/c) 
for some /3n(i) S C°°(IR). The remainder satisfies 

K-1 

(4.19) ||B^(t)||<Cx^||6r^(l,-)||c^(--M+., C;, = CW(M,x,V). 
Hence, we end up with 

v^{l,x) -- 



fe=0 



K-1 



(ZTTfi) 2 \ 



fe=0 



def 1 



-6"(l,a;) + /i^-'^/2i?^(l,x) 



where i?^(l,a:) = (27r)-5 (B^(l, x) +ZYip^„i?^"^(l, •))• Again, 6" satisfies the hypotheses 
of Proposition 10, so we can continue iteratively. To complete the proof, we now have to 
take all the remainders into account. From the discussion above, we get : 

W*P^„v"-i(l,-) = (27rft)-^Zi*P^„6"-i(l,-) + ft^"''/'Z^'P7„(i?r'(l'-)) 

= (27rft)-i (6"(t,-) + S^Bx(i--)) +fi'^"''^'^*P7„(^ir'(lr)) 

where we defined R]^ = {211)- Ib^ + W'P^Ji?^"^). Since |F^,J < 1, we have 

Wu'p^Jl^^l^ <i + ch, 00. 

This implies that ||W'P^„ (i?^-^(l, O)!! < (1 + Cfi)||i?^"i (1, OH, and finally i?^ satisfies 

(4.20) < (1 + ChT {\\Bl\\ + Wk'W +... \\B],\\ + \\Bl\\) 

In view of (4.19) and (4.5), this concludes the proof. □ 

Given w"^^(l, •), we have then constructed v^{t,x) as in (4.15), but it remains to control 
the remainder R^ in norm : from (4.19) and (4.20), we see that it is crucial for this to 
estimate properly the norms of the coefficients for j > 1 and fc € JO, — 1]. 

Lemma 12. Let n > 1, and define 
X>„ = sup 

xe7rA"(l) 



i=0 



X)n = 1. 



If X € 7rA"(l), the principal symbol is given by 

(n-l 
l[v^.,{l,<l>si{x))] 6^(l,0s:(x)). 

For k £lO,K — 1}, the functions b^ satisfy 



26 



EMMANUEL SCHENCK 



(4.22) \mi,-)\\c^ < CkAn + lf+'-Dn 
where Ck,e = C^^'''\M,x,V). It follows that 

(4.23) l|S^(l,-)|| < CKn3^+''l>„_i 

n 

(4.24) mum < CK{i+cnrj2f''+'^-Dj-, 

vjhere C > andCK = C'^^H-^'X;^)- On the other hand, if x ^ 7rA"(l), we have h1{x) = 
for k e 10, if - 11. 

Proof. First, if x ^ 7rA"(l), then there is no p e V^^ such that TT^^{p) = x, and then 
f"(l,a;) = 0{h'^). In what follows, we then consider the case x G 7rA"(l). We first see 
that (4.21) simply follows from (4.18) appHed recursively. If pn = {xn,^n) € A"(l), we call 
Pj = i^i^ij) = *^^"(Pn) e A^(l) if i > 0. In other words, 

Vj G p,n], Xj-i = (l)s^{xj). 

It will be useful to keep in mind the following sequence, which illustrates the backward 
trajectory of p„ € A"(l) under k G |[l,nl and its projection on M : 



eA"-i(l). 



■Pne A"(l) 



Xl 



''Sn 



We denote schematically the Jacobian matrix 



for 1 < I < j < n. Since for 



Sj dxj 

any E > 0, the sphere bundle T*M np~^{E) is transverse to the stable direction [Kli], 
the Lagrangians A" C $"A'' converge exponentially fast to the weak unstable foliation as 
n ^ oo. This implies that $*|ao is asymptotically expanding as f ^ oo, except in the fiow 
direction. Hence, the inverse fiow $~*|a" acting on A" and its projection (/)^* on M have 
a tangent map uniformly bounded with respect to n,t. As a result, the Jacobian matrices 
dxj-i/dxj are uniformly bounded from above : for 1 < i < j < n there exists C = C{M) 
independent of n such that 

dxi 



(4.25) 



dx-i 



< C. 



It follows that if we denote Dj = sup^. T^j{l, xj), there exists C = C{M) > such that 
(4.26) C"^ < Dj < C. 

Note also that 



sup 



[]2?„_,(l,,^-;^(x)) 



0=0 



3=0 



We first establish the following crucial estimate : 

Lemma 13. Let n > 1, and k G |l,n|. For every multi index a of length \a\ > 2, there 
exists a constant Ca > depending on M such that 

d Xn—k 



(4.27) 



dx° 
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Proof. We proceed by induction on fc, from fc = 1 to fc = n. The case fc = 1 is clear. Let us 
assume now that 

d°'Xn-k' 



<Cc,k'"-\ fc'e|l,A;-l] 



and let us show the bound for k' = k. For simplicity, we will denote 

d?"^'^, d-x - 



^ dxf ' ' 5a;«+i ' 

In particular, ||9"a;j|| < Ca- We also recall the Faa di Bruno formula : let 11 be the set of 
partitions of the ensemble {1, |a|}, and for tt e 11, write tt = {Si, ...-Bfc} where -Bj is some 
subset of {1, |a|}. Here |a| > fc > 1, and we denote |7r| = k. For two smooth functions 
g : R'' ^ R'^ and f : R'^ ^ R'^ such that / o g is well defined, one has 

(4.28) dyog^^d^-\f{g)l[d^g. 

Tren Beir 

The term in the right hand side is written schematically, to indicates a sum of derivatives 
of / of order |7r|, times a product of |7r| terms, each of them corresponding to derivatives of 
g of order \B\. It is important for our purpose to note that XI 1-^1 = Continuing from 
theses remarks, we compute 

'= d"Xn-k = dXn-kdnXn-k+1 

+ ^ d^^'^Xn-k Y[ d^Xn-k+l =^ dXn-kXk-l + Yfe-l ■ 
7r6n,|7r|>l Be-TT 

By the induction hypothesis, 

(4.29) < 

since the partitions tt involved in the sum contains at least two elements. Setting M^-i = 
dxn-k , we have 

Xk = Mk-i...MiXi + Mk-2...MiYi+Mk-3...MiY2 

+ ■■■ + MiYk-i . 

From the chain rule we have 

dxj dxj^i+i dxj ' 

and (4.25) yields to \\Mi_i ...Mi\\ = 0{1) for 2 < i < k. Adding up all the terms con- 
tributing to Xk and taking (4.29) into account yields to 

ll^fcll < a(l + 1"-' + 2«-' + . . . + (fc - < Cafc"-^ 

and the lemma is proved. □ 

We now prove (4.22). For this, we will proceed in two steps. First, we show the bounds 
for the principal symbol 6q. Then, we treat the higher order terms b^.k > 1 using the 
bounds on ||&ollc< for a-ny £. For 6q, The C° norm estimate follows directly from (4.21). 
From now on, we denote for convenience 

Vo{xo) = b°{l,xo). 

Computing 

dib^ix^) = diiVr^iXn) . . .V^{xi)Vo{xo)) , 
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we will obtain a sum of terms, each of them of the form 

with a„ + • • • + ao = ^. Note that if £ is fixed with respect to n, most of the multi-indices 
tti vanish when n becomes large : actually, at most \i\ are non-zero, and we will denote 
them by ajj, . . . ,ai^, k < Hence the above expression is made of long strings of T>i , 
alternating with some derivative terms d^'Vi which number depends only on £. We can 
then write 



(4.30) ||M„„...„o||co <fn X 



. . . A, 



Let us examinate each terms d^Vi appearing in the right hand side individually. By the 
Faa di Bruno formula and Lemma 13, we have for i ^ 



(4.31) d^lViixi) = n ^nXi < C^n"-\^\ < 



where = &^^^^\M,xy)- Of course, if i = 0, \\d^Vo{xo)\\co < C„||5^&g||co for some 
constant Ca > 0. Now, for a fixed configuration of derivatives {a} = {ui^, . . -ai^} we have 
to choose ii, . . . ,ih indices among n+1 to form the right hand side in (4.30), and the number 
of such choices is at most of order 0{{n + I)''). Hence, 

Uf^e^nn ^ ^ \- \- \\dT'Vn...dT''-'Vi,_,dT''^iJco 
||<^„&0||CO < T>n 2^ ^ ^ 

{a}H,...,ik 

< D„ 5^ a(n + l)*^(n + l)"i-i . . . (n + 1^"-^ 

{"} 

(4.32) < CfVn{n + lY 

where Ce = C^^^(M,x, V). For higher order terms {b^, k > 0), we remark from (4.17) that 
we can write 

k 

(4.33) bUXn) = V^{XnK-\Xn-l) + E E ^U^n)d^-lK-]i^n-,) ■ 

3 = 1 \a\<2j 

The function F"^ can be expressed with the flow, the damping and the cutoff function F-y„. 
It follows that the norms ||F"^||(7« are uniformly bounded with respect to n: 

l|r7,J|c^ = c(^'^nM,x,v). 

In order to show the bounds (4.22) for fc > 0, we will proceed by induction on the index k. 
The case k = Q has been treated above. Suppose now that for any I and k' € |0, fc — 1] we 
have proven 

\\dibl\\co<Ct{n + lf'''+''D^, Q = C(^'^)(M,x,V). 

As above, to treat the case k' = k, we begin by the situation where £ = 0. To shorten the 
formulffi, we introduce for 1 < i < j < n the functions 

j=l \a\<2j 

JliXn) = Vj{Xj)Vj_i{Xj-i)...Vi{Xi) 
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where the Xi, i < n have to be considered as functions of a;„, namely Xi = 
Iterating (4.33) further, we have : 

bUxn) = J:5ri(x„_i) + r"''= 



n-ibk (a;„_i) + J„r +r ' 



(4.34) = J5'60(.To) + J^Ti''= + j:3T2'^- + ... + j;jr"-i''= + r"''= 

By the induction hypothesis and (4.26), each term r*'*^, i > satisfies 

||r"-^''lico <Cfe(n-^)3'=-lX)„_, 
hence adding up all the terms we get 

ll^^llco < Ck-Dn{bUxo) + ^(n - if-^) < CuVn{n + l)^'^ 

i=0 

and we obtain the bounds (4.22) for £ = 0. To evaluate d^bl, £ > 1, we start from the 
expression (4.34). We notice first that 

/3i+/32=/3 j=l |Q;|<2j 

Using the Faa di Bruno formula and Lemma 13, we get 



^r^-^''(x„-Ollco < C^,i^'-' and ||Sf^5«6^z;-i(x„_,_i)||co < C^j' 



and this implies 

Then, exactly the same strategy used to derive (4.32) shows that 

\\diJ'l+,r'>'\\co < Cen^''-^+^. 
Using these estimates and (4.34) yields to 

WdiblWco < Cein + l)n3'=-i+^ < Q(n + 1)^'=+^ 
where the constant Ce is such that Q = C^^'^\M,x,V). 

□ 

4.5. The main estimate : proof of Proposition 5. As noted before, the Lagrangians 
A" converge exponentially fast as n — ^ oo to the weak unstable foliation. This implies that 
for X € 7rA-'(l), the Jacobians Jsj{x) *== | det (f)g^^^^{x)\ satisfy 

JsAx) 



yj>2, V(x,?)e A^(l), 



1 



< Ce"-''/^, C = C{M) > 0. 



Here, 5" generates the (Lagrangian) local weak instable manifold at point {x,^). Moreover, 
theses Jacobians decay exponentially with j as j — > oo. This means that uniformly with 
respect to n, 

n—l n—1 

n Js„-A&^)) < cm n Js-(*-.- (.^sf (a;)) . 

The Jacobian Js^{x,£,){x) measures the contraction of along the unstable subspace 
i;"($i(p)), where $i(p) = (a;,^), and a; e M serves as coordinates to compute this Jacobian 
(via the projection tt). The unstable Jacobian J"(p) '= | det (rf^~^|B",o($(p))) | defined in 
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Section 2.1 express also this contraction, but in different coordinates: for n large enougfi, 
tfie above inequality can then be extended to 

n — 1 n— 1 n — 1 

(4.35) n Js^.,{^s'M) <C l[ Js.(<i>-. (</>s;:(x)) < C [] J"(<i>-^(p)) . 

where C, C only depends on M. As noted above, because of the Anosov property of the 
geodesic flow, the above products decay exponentially with n. Together with the fact that 
the damping function is positive, it follows that the right hand side in (4.23) also decay 
exponentiahy with n. Recall now that 1 < n < Nto and N = Tlogh~^. Using (4.24), we 
then see that the remainders in (4.15) are uniformly bounded : they satisfy 

\m\<CK, Ck = CW(M,x,V) 

uniformly in n and zq, the point on which was based. From the very construction of 
h^{t,x), we then have 

(4.36) \\U'P^„ . ..U'P^^UH^ - ^^-1^6"(1, Oil < Ck/i^-'^/'. 

But the bounds on the symbols 6^, fc > given in Lemma 12 tehs us that (4.36) also holds if 
we replace the full symbol 6" by the principal symbol 6q , provided h is chosen small enough 
- say h < ho{e). Hence, for h < ho, 

\\l{^P^^...U^P-y,l{^6J < {2TTh)-i\\b^{l,-)\\+CKh^-'^^^. 
Now, using (4.21), (4.35) and the fact that |F^| < 1, we conclude that for a" as in (1.8), 

n 

||6^(l,x)|| < Ce"^(^) sup exp^a"o$-^(x,d,5„(l,:r)) 

Here, C = C{M) depends only on the manifold M. Let us consider now the particular case 
n = Nto with N = T\ogh~^. It fohows immediately that 

Nto N / to-1 

sup exp a" o ^~^{x, dxSn{l, x) <Y\ sup cxp a" o $-' (p) 

By the superposition principle already mentionned in (4.4), we then obtain for some C 
C{M) > depending only on M: 

\\UP^,,^...UP^M^O^f,{x)\\ < Cj2sup\\UP^,,^^...P,,U%J\ 

N / io-1 

< Ch-iJ] sup expVa"o$J(p) 

To get the last lino, wo have noticed that K can be chosen arbitrary large: since n < 
Tto^ogh~^ , we see that for h small enough, the main term in the right hand side of the 
second line is larger than the remainder CKfi^~'^^'^ , and e-'^*"'-'''') = 0{1). This completes 
the proof of Proposition 5. 
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Appendix A. Semiclassical analysis on compact manifolds 

In this appendix we gather standard notions of pseudodifferential calculus on a compact, 
d dimensional manifold M endowed with a Riemannian structure coming from a metric g. 
As usual, M is equiped with an atlas {fe, Ve}, where {Vt} is an open cover of M and each 
fe is a diffeomorphism form Ve to a bounded open set We C M**. Functions on R'^ can be 
pulled back via /; : C°°{We) C°°{Ve). The canonical lift of fe between T*Ve and T*We 
is denoted by ff. 

{x,0 G T*Ve ^ fe{x,0 - {fe{x),{Dfe{x)-YO e T*We , 

where denotes the transpose of A. Its corresponding pull-back will be denoted by 
/; : C°°{T*We) C°°(T*V«). A smooth partition of unity adapted to the cover {Ve} is a 
set of functions 4>e € C^{Ve) such that ipe = 1 on M . 

Any observable (i.e. a function a G C°°(T*M)) can now be split into a = J2e'^^ where 
ae = (pea, and each term pushed to ag = G C°°{T*We). If a belongs to a standard 

class of symbols, for instance 

a e S"'''' = S'^iiO"') =^ {a = Oft e C°°(M), \d^d^a\ < ^^^^/^-'^(O'"-""} , 

each ae can be be Weyl-quantized into a pseudodifferential operator on S(R) via the formula 

Vu e S{R% Op^{de)u{x) = I e*<--^'«> ae ^; u{y) dy d^ 

To pull-back this operator on C°°{Ve), one first takes another smooth cutoff ipe € C^{Ve) 
such that = 1 in a neighbourhood of supp 4>e . The quantization of a e is finally 

defined by gluing local quantizations together, yielding to 

V« e C-(M),Opft(a)u = X /; o Op^(a,) o {f^'n^eu) 

I 

The space of pseudodifferential operators obtained from 5''=''" by this quantization will be 
denoted by v]/™'*=. Although this quantization depends on the cutoffs, the principal symbol 
map a : gmM j gm,,k-i jg intrinsically defined and do not depend on the choice of 

coordinates. The residual class is made of operators in the space \1>™'~°°. As an example, the 
(semiclassical) Laplacian —h^Ag e ^""'^ is a pseudodifferential operator, and its principal 
symbol is given by a{-h^Ag) = = g,(^,^) e S^^°. 

In this article, we are concerned with a purely semiclassical theory and then deal only with 
compact subsets of T*M. If A G we will denote by WFft(^) the semiclassical wave 

front set of A. A point p e T*M belongs to WF ti{A) if for some choice of local coordinates 
near the projection of p, the full symbol of A is in the class WFft(yl) is a closed 

subset of T*M, and WF?i(AB) C WFft(A) n WFft(B). In particular, if WFft(A) = 0, then 
j4 is a negligible operator, i.e. A € If \1> g L'^[M), we also define the semiclassical 

wave front set of ^ by : 

WFfiX*) = {{x,i) ■■ 3a e a{x,0 ^ 0, || Op,(a)vI/|U2(M) = 0{h^)Y 

where the superscript ° indicates the complementary set. We will often make use of the 
following fundamental propagation property : if Z7* is a Fourier integral operator associated 
to a symplectic diffeomorphism $* : T*M T*M, then 



WFft(f/**) = $*(WFft(*)). 
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